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. Abstract. The Landau Hamiltonian, describing the behavior of 

a quantum particle in dimension 2 in a constant magnetic field, is 
perturbed by a magnetic field with power- like decay at infinity and 
CO ' a similar electric potential. We describe how the spectral subspaces 

change and how the Landau levels split under this perturbation. 
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1. Introduction 



The paper is devoted to the study of the spectrum of the Schrodinger 
and Pauh operators in the plane, with nonzero constant magnetic field 
(N ■ perturbed by a smooth magnetic and electric fields that decay power- 

like at infinity. It continues the paper [15] by the authors, where the 
^ ', case of a compactly supported perturbation was considered. 

O I The Landau Hamiltonian describing the motion of the quantum par- 

ticle in two dimensions under the influence of the constant magnetic 
field is one of the classical models in quantum physics. The spectrum 
0\ ', (found first in [4], see also [9]) consists of eigenvalues with infinite 

^ I multiplicity lying at the points of an arithmetic progression. These 

eigenvalues are traditionally called Landau levels (LL) and the corre- 
^ . sponding spectral subspaces are called Landau subspaces. 

^ I A natural question arises, what happens with the spectrum of the 

Landau Hamiltonian under the perturbation by a weak electrostatic 
potential or/and magnetic field. One should expect that the Landau 
levels split, and the problem consists in describing quantitatively this 
splitting as well as in studying the behavior of the spectral subspaces 
under the perturbation. 

The case of the perturbation by an electric potential V, moderately, 
power-like decaying at infinity, was first studied by Raikov in [11]. The 
case of a fast decaying (or compactly supported) electric potential was 
dealt with much later, in [13], see also [10]. It was found that the 
character of the Landau levels splitting depends essentially on the rate 
of decay of V and the asymptotics of the eigenvalues in clusters can be 
expressed in the terms of V, quasi-classically or not. 
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When the magnetic field is also perturbed, the situation becomes 
more complicated since is not the magnetic field itself but its potential 
that enters in the quantum Hamiltonian. So, the perturbation of the 
operator turns out to be fairly strong even for a compactly supported 
perturbation of the field and it may even be not relatively compact if 
the perturbation goes to zero at infinity not sufficiently fast. Iwatsuka 
[6] proved that the invariance of the essential spectrum still takes place, 
so Landau levels are the only possible limit points of eigenvalues lying 
in the gaps between them. Further on, in [3], [7], [8], [12] the char- 
acter of the splitting of the lowest Landau level was investigated. For 
compactly supported magnetic field perturbation and electric potential 
the splitting of all Landau levels was studied in [15]. The main result 
of [15] was the description of the spectral subspaces of the perturbed 
operator corresponding to the clusters around Landau levels. It was 
found that these subspaces change fairly strongly, and a rather exact 
approximation for these subspaces was found in the terms of modi- 
fied creation and annihilation operators. At the same time, although 
the perturbation of the operator may be very strong, the splitting of 
eigenvalues is super-exponentially weak, just like it is in the case of a 
perturbation by a compactly supported electric field only. 

In the present paper we continue the study of the Landau Hamil- 
tonian with a perturbed magnetic field. Now we consider the case of 
the perturbation decaying moderately, power-like, at infinity. For the 
spectral subspaces the results similar to the ones in [15] hold. As for 
the asymptotics of the eigenvalues in clusters, we obtain more com- 
plete results, proving estimates and, under some natural conditions, 
the power-like asymptotic formulas for eigenvalues. The problem of 
the splitting of Landau levels under moderately decaying perturba- 
tions of the magnetic field has been considered in [5], see Theorem 
11.3.17 there, where even the remainder term of the asymptotics of 
the eigenvalues in clusters was found. Powerful methods of microlocal 
analysis are used in [5]. Our methods are more elementary, give the 
approximation for the spectral subspaces, moreover, we do not need 
the rather restrictive 'hyperbolicity' condition (11.3.49) imposed in [5]. 

The paper is heavily based upon the methods and results of the 
papers [7] and [15]. Following the structure of the latter paper, we 
refer rather shortly to the fragments that should be repeated, word for 
word or with minor changes only, in order to be able to concentrate on 
important differences. For convenience of references, the numbering of 
sections here coincides with that in [15]. 

The second author was partly supported by a grant from the Royal 
Swedish Academy of Sciences, he also thanks Chalmers University of 
Technology in Gothenburg for hospitahty. 
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2. Magnetic Schrodinger and Pauli operators 

2.1. The unperturbed operators. We will denote the points in the 
plane by x = {xi,X2); it is convenient to identify with C by 
setting z — xi + 1x2- So, the Hilbert space L2{^) with Lebesgue 
measure (which will be denoted by dx) is identified with L2(C). The 
derivatives are denoted by by dj- — d^^. and we set, as usual, d — 

{di + 182) /2, d = {di- 182) /2. 

The constant magnetic field is denoted by B° > 0. The correspond- 
ing magnetic potential is A'{x) — (74^,^42) = ^{—X2,xi). Then the 
(unperturbed) magnetic Schrodinger operator in L2(]R^) is 

H° = -(V + a°)^. (2.1) 

The Pauli operator describing the motion of a spin-| particle acts in 
the space of two-component vector functions, it has the diagonal form, 
P° = diag (P;, P° ), where P;^ = H° ± B°. 

The spectrum of the Schrodinger operator is described by the classi- 
cal construction originating in [4] . For the complex magnetic potential 
A° — A\-\- 2^2 the creation and annihilation operators are introduced, 

= -2id -A^, Q° = ~2id - A°. (2.2) 

These operators can be also expressed by means of the scalar potential, 
the function "^°[z) = solving the equation A\E'° = B°: 

Qo = -2ze-*°ae*°, = -2ze*°ae-*°. 
The operators Q°, Q° satisfy the following basic relations 

[Q°,Q^] = 2B°. (2.3) 

P^ = Q°Q^, P^ = Q^Q°, H° = Q°Q^ - B° = Q^Q° + B°. (2.4) 

The spectrum of H° is described in the following way. The equation 
P'Lu = 0, M G L2 is equivalent to Q°m = e^*°c?(e*°M) = 0. 

This means that / — e^°u is an entire analytical function in C, such 
that after being multiplied by it belongs to L2. The space of such 
functions / is called Fock or Scgal-Bargmann space JF. 

So, the null subspace of the operator PI, i.e., its spectral subspace 
corresponding to the eigenvalue Aq = 0, is Cq = e~^° T . After this, 
by the commutation relations (2.3), (2.4), Lq — Q°>Co are the spectral 
subspace of PI with eigenvalues Ag = 2gB°, g = 0, 1, . . . , called Landau 
levels, and the spectra of H°, P^|_ consist, respectively, of A^ + B° 
and Ag + 2B°. The operators Q°,Q° act between Landau subspaces 
= Q^^£o,g = 0,l,..., 

W-J^g^ Q° : A ^ Q° : >Co ^ 0, (2.5) 

and are, up to constant factors, isometries of Landau subspaces. 
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2.2. The perturbed operator. We introduce the convenient class of 
functions. A function F{x) is said to belong to the class Sjj, /3 < 0, if 



F(x) = 0{\xf), dlF{x) = 0{\x\P-'), oo, = 1, 2, . . . , (2.7) 



for some 5 e (0,-/5). The particular value of 5 is irrelevant. So, if 
F e Sp then djF can be considered as a function in Sp_s (with | 
acting as S). 

Now wc introduce the perturbation b G C°°(]R^) of the magnetic 
field and set B = B° + b. We suppose that 



Under the condition (2.8), the upper estimate for the counting func- 
tion of the eigenvalues in the clusters will be proved, as well as the 
approximate representation of the spectral subspaces, following [15]. 
For obtaining asymptotic formulas, additional conditions will be im- 
posed on b and eventually on the electric perturbation. 

Let ■0 be a scalar potential for the field b, a solution of the equation 
Ai/j = h. Of course, ip is defined up to a harmonic summand, the 
choice of ip corresponds to the choice of gauge. The magnetic potential 
a = (01,02) = {—d2ip,diip), curia = b, is thus determined up to a 
gradient, and the complete scalar and vector magnetic potentials are 



We define the perturbed magnetic Schrodinger operator as in (2.1), 
with ^4° replaced hy A: H = —{V + iA)^, and the components of Pauli 
operator as P± = H ± B. It is easy to observe that the difference 
between H and H° contains an operator of multiplication by A° ■ a, the 
latter function does not decay at infinity, and thus looks like being not a 
relatively compact perturbation of H°. However, thanks to the special 
form of this term, the perturbation is still relatively compact if a at 
infinity, as was noticed by Besch [1]. In our case, under the conditions 
imposed on b, the scalar potential grows at most logarithmically, and 
the vector potential a decays I at infinity. 

Now let us look at the algebraic structure related with the perturbed 
operators. The perturbed creation and annihilation operators are de- 
fined similarly to (2.2): Q = —2id — A, Q = —2id — A, where A is the 
complex magnetic potential, A = A° + (oi -|- ^02) . The commutation 
relations for Q, Q have the form 



h e Sfs for some /3 < —2. 



(2.8) 



^ = ^° + ^, A = A° + a. 



(2.9) 



[Q, Q] = 2B = 2B° + 2b, 



(2.10) 
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and P± and H satisfy 

P+ = QQ, P_ = QQ, P+ - P = 2B = 2B° + 2b, (2.11) 
H = QQ-B = QQ + B. (2.12) 

The relations (2.10)-(2.12) contain variable functions on the right- 
hand side and the spectra of Schrodinger and Pauli operators do not 
determine each other any more. The only information that one can 
obtain immediately, is the description of the lowest point of the spec- 
trum of P_. Since, again, Q = Q*, the equation P_m = is equivalent 
to Qm = 0, or (9(exp(\E')-u) = 0. So the function / = -uexp^I/ is an 
entire analytical function such that u = exp(— ^)/ e L2. The space of 
entire functions with this property is, obviously, infinite-dimensional, 
it contains at least all polynomials in z variable, although it does not 
necessarily coincide with the Fock space. We denote the null-space of 
P_, the space of zero modes, by Tio. It is infinite-dimensional; complex 
polynomials times exp(— \l/) form a dense set in Hq. The lowest Landau 
level Ao is an isolated point in the spectrum of P_. 

As it follows from the relative compactness of the perturbation, by 
Weyl's theorem, the essential spectrum of the perturbed operator P_ 
consists of the same Landau levels Ag, and the eigenvalues in the gaps 
may only have A^ as their limit points. This latter fact was estab- 
lished much earlier by Iwatsuka [6], and the constructions in [15] can 
be considered as the extension of the approach in [6]. 

Now we add a perturbation by the electric potential. Let V{x) be a 
real valued function in 5/3, /3 < —2. We introduce the operators 

H{V) = H + V, P±{V) =Pi, + V. (2.13) 

Since the operator of multiplication by V is relatively compact with 
respect to H, P±, the operators (2.13) have the same essential spectra 
as the respective unperturbed ones (Aq ceases to be an isolated point 
of the spectrum of P-{V) ). In this paper we are going to study the 
distribution of the eigenvalues of il{V), P±(^) near A^. 

2.3. Some resolvent and commutator estimates. In the course of 
our proof we will need some boundcdness property of the resolvent of 
the operators P+,P_,H and their spectral projections. These prop- 
erties, in a slightly less general form, were established in [7], Lemma 
1.4. 

Wc denote by R±{z) the resolvent of the operators P± and by Ilj, 
j = 1, 2, the operators 11^ = idj + {Aj + a^). 

Proposition 2.1. Suppose that V,h belong to S^, u < —2. Let P he 
the spectral projection of the operator P± or H corresponding to some 
hounded isolated piece of the spectrum. Then for any real I the fol- 
lowing operators are hounded: (x)~^+'^~'[y, P](x)', (x)~^"'"'^~'[b, P](x)', 
{x)-^R±{xY, {x)-%R±{xy. 
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The proof in [7] is given for the class with 5=1, for the projection 
P corresponding to the lowest Landau level and for positive I only. In 
our formulation, the proposition is proved in an analogous way. 

Next wc establish the estimates for eigenvalues and singular numbers 
of some compact operators. The estimates will be needed further on, 
in the process of proving the required eigenvalue asymptotics. 

For a compact operator T we, as usual, denote by n(A, T) the dis- 
tribution function of the singular numbers (s-numbcrs) of T, i.e. the 
quantity of s- numbers of T that are bigger than A. If the operator 
is self-adjoint, the distribution functions for the positive and negative 
eigenvalues of T are denoted by n±(A, T). The operator T can be 
dropped from the notation if this does not cause misunderstanding. 

Proposition 2.2. Let V he a function in Si,, v < Consider the 
operator X — X^iy) — VP^^ . Then for N sufficiently large, 

n{X,X)^0(X^), A^O. (2.14) 
Moreover, if K is a compact operator then 

n{X,KX), n{X,XK) = o{X^), A ^ 0. (2.15) 

Compared with Proposition 5.1 below, the above estimate shows 
that for N large enough, the operator Xn{V) admits the same spectral 
estimates as the Toephtz type operator, with P+'^ replaced by the 
spectral projection of P±. 

Proof. Consider a very special case first. Let V e Si,,u e (— |, 0) and 
N e (— z/, |]. We will obtain the s-numbers estimates of the oper- 
ator VP^^ . These numbers are majorated by the s-numbcrs of the 
analogous operator with V replaced by = (x)"^'^. The semigroup 
generated by P_— iV, by the diamagnetic inequality, is dominated by 
the semigroup generated by (1 — A) ^, therefore, by the results of [14] 
(see Theorems 1, 4 and Sect. 5.4 there), these singular numbers are ma- 
jorated by the s- numbers of the operator Tjy = Vp{l — A)"". By the 
Fourier transform, this operator is unitary equivalent to (1 — A)''(^)~^, 
and for the latter operator the required eigenvalue estimate is given by 
Cwikel's theorem in [2]. 

For the general case, again lor V — Vj,, v < 0, we take N so large 
that -u/N < \. We represent V in the form V = V = V^+i. 

Now in the expression \/^+ip^^ -^ve leave one copy of V in the first 
(utmost left) position and start moving the remaining copies to the 
right, commuting them with copies of P^^ in such way that finally 
there will be only one entry of V or its derivatives between two copies 

of P+^. As a result, we arrive to a collection of summands, each being 

_i _i 

the product of operators of the form l^P^^, P^^l^, and, possibly, 
some more bounded operators of the form considered in Proposition 
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2.1, where W G S'^/at is V or some of its derivatives. For the operators 
_ 1 _ 1 

WP_^_^, P_|_^iy we can apply the estimate found in the first part of 
the proof and obtain the required inequality using the Weyl inequality 
for the s-numbers of the product of operators. The second statement 
follows from the first one and, again, the Weyl inequality. □ 

3. Approximate spectral subspaces 

In this section, under the condition that b satisfies (2.8), we con- 
struct the approximate spectral subspaces of the operators (2.13). This 
is done in the same way as in [15], so we just briefly describe the con- 
struction and pinpoint the main differences. 

First of all, we consider the null subspace Tio. In [15] it is shown 
that for b G C^, Hq possesses a dense subspace of rapidly decaying 
functions. The same reasoning proves this property for our case. It is 
here, that the condition /3 < — 2 implying that V, b G Li is essential. 

Let = (Aq — 7, Aq + 7), g = 0, 1, 2, . . . , 7 < B°, be intervals of the 
same size on the real axis, centered at the Landau levels A^ = 2gB°. 
We choose the size of 6q in such way that neither of these intervals has 
the eigenvalues of P_ at its endpoints. Moreover, since the lowest LL 
Ao = is an isolated point of the spectrum of P_ , we can choose the size 
of the intervals in such way that So contains only this point of spectrum. 
We denote by Tig the spectral subspace of P_ corresponding to the 
interval 6g and by Pg the corresponding spectral projection. Since, by 
[6], the spectrum of P_ is discrete between Landau levels, the change of 
5q leads only to a finite-rank perturbation of Pg. As usual, the spectral 
projection Pg can be be expressed by means of the integration of the 
resolvent of P_ along a closed contour Tg in the complex plane, not 
passing through the eigenvalues of P_ and containing inside only those 
eigenvalues that lie in 5q. Again, using the discreteness of the spectrum 
of P_ between the Landau levels, we can choose these contours so that 
they are obtained from Fq by the shift along the real axis in the complex 
plane, F^ = Fq + 2qB°. 

Now we are going to estabhsh several properties of the subspaces Hg, 
projections Pg and some related operators. First, note the simple fact 
following directly from the spectral theorem. 

Proposition 3.1. For any q — 0,1, ... , and any polynomial p(A) the 
operator p{P ^)Pg is bounded, moreover (p(P_) —p{Ag))Pg is compact. 

In fact, by the spectral theorem the nonzero spectrum of the operator 
p(P_)Pq consists of the points p{Xj) where Xj are all points of spectrum 
of P_ in 6g and thus all p{\j) live in a bounded interval. Moreover, 
p{\j) may only have p{Ag) as their limit point by Iwatsuka's theorem. 

The following lemma will enable us later to prove a much stronger 
compactness property. 
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Lemma 3.2. Let each ofTj, j = 1, . . . ,N be one of operators Q or Q. 
Then for some constants C, C , for any u in the domain of the operator 

-pN 

llTiTs . . . TnuW^ < C(P%, u) + C'\\u\\\ (3.1) 

The proof of Lemma 3.2 for the case b e can be found in [15] 

Section 7, in our case the proof goes exactly in the same way. 
Now we can estabhsh the compactness property. 

Proposition 3.3. Let Ti, . . . ,Tn be a collection of operators, each be- 
ing Q or Q, and let hj , j = 0, . . . N be functions with all derivatives 
bounded, T — hoTihi . . . Tjsihiq. Then for any q and for any polynomial 
p{X) the operators T(p(P_) — p{Aq))PqT are compact. 

Proof By commuting functions hj and operators Tj (moving all func- 
tions to the left), we transform the left operator T to the sum of terms 
of the form H^Tk where are bounded functions and 7^ is a product 
of no more than N operators Q, Q. Similarly, in T that stands to the 
right of Pg, we move all functions to the utmost right positions, to get 
the representation of T as a sum of terms T^^h^, being the product of 
a bounded smooth function h:^ and no more than N operators Q, Q. 

For each of the terms arising in this way in the decomposition of 
T(p(P_) — p{Ag))PqT, we can write 

KTMP-) -p{K))PgTji^ = [ha.{^- + l)-']x (3.2) 

[(P^ + l)(p(P_) -p(A,))P,(P^ + l)][(P!f + 1)-'TM- 

In (3.2), the first factor in brackets is bounded by Lemma 3.2, and the 
middle factor is compact by Proposition 3.1. The last factor in brackets 
is also bounded, by Lemma 3.2 applied to the adjoint operator. □ 

Now we describe the main construction of the paper, the approximate 
spectral subspaces of the perturbed operator. It is sufficient to consider 
the operator V^{V). In fact, by (2.10), (2.11), (2.12), 

H(\/) = P_(\/ + b) + B°, P+(\/) = P_(\/ + 2b) + 2B°, (3.3) 

and thus these operators differ from P_(V^) by a shift and by the electric 
type perturbations b, 2b. We find the approximate spectral subspaces 
of P_. Adding an electric perturbation will then be an easier task. 
The subspaces approximating Tig will be defined as: 

Q^^Uo, Q.^Q'Go: g = l,2,.... (3.4) 

So we mimic the construction of the eigenspaces of the unperturbed 
Landau Hamiltonian, see (2.5), in the same way as it was done in [15], 
by applying the creation operators to the space of zero modes. 

Of course, since we apply the unbounded operator Q, we must show 
that we never leave the space L2, and moreover, that the subspaces Qg 
are closed. Both these properties, as well as some other results will be 
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based upon the important Proposition 3.4 (an analogy of Proposition 
3.4 in [15]). 

The essential difference is that now, when the perturbation docs not 
have compact support, wc have to trace the rate of decay of different 
terms arising in the process of transformations. This analysis enables 
us to single out the leading terms in the resulting expansions. 

Proposition 3.4. Let q > 0. 

(1) There exists a function Zq[h] G Sp depending only on q, B°, 
and b such that for any u & Tto, 

WQ^uf = Cq\\u\\^ + {Zg[h]u,u), Cg = q\{2By. (3.5) 

The function Zq[h] is a polynomial in b and its derivatives up 
to the order 2q — 2 with coefficients depending on B°. The term 
linear in b and not containing derivatives equals C'gB°'^'~^h, 
Cg = 2«g!g. Moreover, 

Zq[h] - C^B°«-^b = 0{\xf-^) (3.6) 

at infinity. 

(2) Let U{x) be a function in Sp,(5 < —2. There exists a function 
Xq\b, U] e Sp depending only on q, B°, b, and U such that for 
any u e TYo- 

(C/QV Q%) = (A'Jb, U]u,u). (3.7) 

The function A'^fb, C/] is expressible as an order 2q linear dif- 
ferential operator acting on U, with coefficients depending poly- 
nomially on b, its derivatives, and B°, moreover, 

Xg[h, U] - C'gB°m = 0(1x1''-^) (3.8) 

at infinity. 

Proof. The combinatorial part of the proof, consisting of multiple com- 
muting of the creation and annihilation operators with functions and 
with each other is exactly the same as in [15]. What remains to be 
checked are the estimates (3.6) and (3.8). These estimates follow from 
the fact that all terms in the expressions -Zg[b], A'g[b, ?7], except the 
leading ones, contain either derivatives of b or U, or products of these 
functions, and therefore decay at infinity not slower than □ 

When applying Proposition 3.4 and similar results, we need a certain 
compactness property. Such facts were used persistently in [13], [10], 
but for the case of a constant magnetic field only. 

Lemma 3.5. Let W{x) be a function in S^, < 0. Let £ be an 
arbitrary differential operator having the form 



SL^flTj2T2...Trn, 



(3.9) 
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where each of Tj is one of operators Q, Q and fj are functions, with 
all derivatives bounded. Then the quadratic form 

to[u] = J W{x)\£ufdx (3.10) 

is compact in the space Ho- 

Proof. Let u — e~*/i be a function in Ho, so h{z) is an analytical 
function. We write the quadratic form (3.10) as 

to[u] = {W£u,£u). (3.11) 

Now we move all the operators Q in £ from the second factor in (3.11) 
to the first one and from the first factor to the second one; thus they 
turn into Q. In the process of commuting these Q with Q and the func- 
tions W and fj, some derivatives of these functions appear; the function 
W goes to zero at infinity, together with derivatives, the derivatives of 
fj are bounded. Then, by means of commuting the operators Q with 
the functions, we move all entries of Q in the first and in the second fac- 
tors in (3.11) to utmost right position, where they vanish since Qu — 
for u e Tio. The only remaining term in the form tt3[u] will be 

tt)[u]^ {Wiu,u), (3.12) 

where Wi is a function (a combination of W.fj.h and their deriva- 
tives) tending to zero at infinity. Now take e > and represent Wi as 
Wi = Wi^s + W^i',£ so that \W[^J\ < e and Wi^^ has compact support. 
For (W[^^u,u), we have the estimate by so the corresponding 

operator has norm not greater than e. For (l^i^it, it), we take some 
R such that the support of the function Wi^^ lies inside the circle Cr 
with radius R centered in the origin. For each r G (R, 2R) we write 
the Cauchy representation for an analytical function h{z): 

h{z) = {2mY' h{C){z - Cr'dC. (3.13) 

for some fixed function ^(r) G C^(R, 2R), J ^{r)dr = 1, we mul- 
tiply (3.13) by ^{r) and integrate in r from R to 2R. This gives 
the integral representation of h{x), \x\ < R, in the form h{x) = 
J-^^^^^^^^K{x,y)h{y)dy, with smooth bounded kernel K{x,y). Af- 
ter applying £ in x variable, we obtain the representation for £u — 
£(e-*/i): 

2u{x) = / e-*(^)X^(a;, y)e''^y^u{y)dy = {K.^u){x). 

-'R<|2/|<2R 

The integral operator iH^i^gl^/^K"^ has a bounded kernel with compact 
support and therefore is compact in L2, and thus the quadratic form 
Vo[u] can be written as ro[ii] = (sign VFi^ejl^i^el^/^K^w, |1^i^£|^/^K^m) 
and therefore it is compact. Now we see that the quadratic form (3.10) 
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can be for any e represented as the sum of a compact form and a form 
with norm less than £, and this proves the required compactness. □ 

Now we are able to justify our construction of the spaces Qq. 

Proposition 3.6. The sets Qq defined in (3.4) are closed subspaces in 

Proof. The fact that Qq C L2 follows directly from Proposition 3.4. 
Next, the relation (3.5) can be written as 

(PoQ'^Q''«,m) = Cq{u,u) + {PQZq[h]u,u)- u = Pm e Ho. (3.14) 

In the second term in (3.14), by Lemma 3.5, the operator Po-2g[b] is 
compact in Hq, and therefore we can understand (3.14) as showing 
that the operator C~^PoQ^ is a left parametrix for Q"^ : Ho ^ L2. 
This imphes that the range of is closed. □ 

The mill space of Q and therefore of Q'' is zero. Consider the op- 
erator Q*^ as acting from Qq = TYq to Qg. This is a bounded invertible 
operator, therefore the inverse, that we denote by Q~^, is a bounded 
operator from Qg to Qq. It is a compact perturbation of Q^. 

4. Approximate spectral projections 

In this section we prove that the subspaces Qq are very good ap- 
proximations to the spectral subspaces Tig of the operator P_, and to 
the spectral subspaces of P_(V). Closeness of subspaces will be mea- 
sured by closeness of orthogonal projections onto them. Recall that 
the projection onto Hq is denoted by Pq. Let Qq be the projection onto 

Theorem 4.1. The projections Pg and Qq are close.' for any N, and 

any collection of the operators Tj, j = 1, . . . ,N, each ofTj being Q or 
Q, the operator T{Pq — Qq)T , is compact, T = T1T2 . . . T^q. 

In justifying the theorem, we need two technical lemmas, both con- 
cerning the properties of products of many copies of the resolvents of 
P+ and P_, the creation and annihilation operators Q and Q, functions 
hj with all derivatives bounded, and, possibly, the spectral projection 

Pq. In such product, we assign order 1 to Q and Q, order —2 to the re- 
solvent, order to functions and projections. The order of the product 
is defined as the sum of orders of factors. 

Lemma 4.2. Let 21 be the product of creation, annihilation operators, 
resolvents, and functions hj, have negative order, and let at least one 
of the functions hj belong to S^,, ly <0. Then 21 is compact. 

Lemma 4.3. Let 21 be the product of creation, annihilation operators, 
resolvents, functions and the projection Pg. Then 21 is bounded. If, 
moreover, at least one of hj belongs to Si,, u <0, then 21 is compact. 
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The proof of the earher versions of these lemmas, with the Sy - 
condition replaced by the compactness of the support of one of the 
function is given in [15]. In the present formulation, the proof goes 
exactly in the same way, only using the new version of Lemma 3.2 and 
Lemma 3.1, fit for the relaxed conditions for the functions. 

The proof of the theorem goes in the following way. We construct an 
intermediate operator &q with range in Qq and prove that &q is close 
both to Pq and Qq. For > we define the operator &q as 

&q = C-^Q^PoQ^ Cq = q\{2^y. (4.1) 

So, our expression for the approximate spectral projection is just a 
natural modification of the exact formula (2.6) for the unperturbed 
operator. Equivalently, the operator <Sq can be described by the for- 
mula 

e, = C-^GG*, G = Q^Po; 

in Proposition 3.4 this operator is shown to be bounded. 

The proof of Theorem 4.1 will consist of two parts, showing that &q 
is close to Pq and showing that it is close to Qq. The second part is 
proved exactly like in [15], since it is based upon Lemma 4.3 only. As 
for the first part, we need a more detailed information of the difference 
&q — Pq. This information is given in the following statement. 

Proposition 4.4. The operator &q is close to the projection Pq. More- 
over, the difference &q — Pq has the the form 

&<i-P,- C.CP- - ^,)P, + Z„ (4.2) 

where Zg is such an operator that {x)^^^^TZqT" is bounded for any 
T, T' being finite products of creation and annihilation operators. 

So, the improvement, compared with Proposition 4.4 in [15], consists, 
first, in the separation of the leading term in the difference, the operator 
(P_ —Aq)Pq, and, secondly, in the establishing the improved smallness 
of the remainder term Zg 

Proof. Combinatorically, the reasoning goes in the same way as in [15] , 
with a natural replacement of the auxiliary results requiring the com- 
pactness of the support of the perturbation and singling out the leading 
term. We remind the main structure of the reasoning, omitting the de- 
tails. The case g = 1 is considered first; this case contains all typical 
features. The higher Landau levels, q > 1, can be taken care of in the 
same way as in [15], by the induction on q. 

Recall that R±{C) denotes the resolvent of the operator P-t. The 
projection Pq can be expressed via Riesz integral 

Po = {2mr' [ R-iOdC, 

J To 
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where Tq is the closed curve defined, together with curves Tq, in Sect. 3. 
We are going to transform the expression for the resolvent -R- (C) using 
the commutation relations (2.10), (2.11), (2.12). After this, the crucial 
observation is that the integral of i?^, A- > 2 along Tg vanishes. This 
enables us to dispose of terms that before integration were not weak 
enough. We start by writing 

i?_(C) = (P+ - 2B° - 2b - C)-' = (P+ - 2B° - C)"'- (4.3) 

(P_ - C)"'(2b)(P+ - 2B° - C)-' = i?+(2B° + C) - Z{C) 

We multiply (4.3) by Q from the left and by Q from the right, as (4.1) 
requires. For the first term we use that 

Qi?+(2B° + C)Q = Q(QQ - 2B° - C)"'Q = (4.4) 

QQ(QQ - 2B° - C)"' = P-R-{2B° + ()■ 

Integration gives 

/ Q(P_ - 2B° - C)-'QrfC = P- / (P--C)"'^^C=(27ri)P-Pi. 

So, we have 

Si - Pi = A-^(P - Ai)Pi + Zi, (4.5) 

where Zi is the contour integral of the term Z{() in (4.3). 

As explained in Section 2, the nonzero eigenvalues of P_Pi may con- 
verge only to Ai, and thus (P_ — Ai)Pi is compact by Proposition 3.3. 
Moreover, this compactness is preserved after the multiplication by any 
product of creation and annihilation operators. 

Now we consider the remainder term, Z(() in (4.3). We are go- 
ing to show that for any operators of the type T, T' the integral of 
{x)~^^^T Z{QT' along Fq is bounded. Due to the arbitrariness of T, T' 
this , of course, implies the compactness required by the theorem. 

Let be some fixed, sufficiently large integer. We apply the resol- 
vent formula (4.3) 2A^ times to the first factor in Z[C,). This operation 
will produce terms of order —4, —6, . . . , containing factors R-{C) and 
b, and the remainder term of order — 2A^ containing these factors with, 
additionally, one factor P+(2B° -|- (). This last remainder term, for A^ 
sufficiently large, is compact by Lemma 4.2, even after the multiplica- 
tion by T, T'. We will study its spectral properties below. 

The leading terms in Z{(), having orders —4, —6,. . . , will be trans- 
formed by repeatedly commuting b and -R-(C) and then the resulting 
commutants again with R-{C) and so on. Under commuting R-{C) 
with a function, with Q, or with Q, this factor P-(C) rnoves to the left 
or to the right, and one more product in the sum composing TZ{C,)T' 
arises, of the order lower by 1. In this commuting procedure we aim 
for collecting the factors P_ (C) together. As soon as we obtain a term 
with all P_ (C) standing together, we leave it alone and do not trans- 
form any more. After sufficiently many commutations, we arrive at a 
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collection of terms of negative order, smaller than —N, in TH{QT'. 
all of which are compact by Lemma 4.2. The terms of order —N or 
higher will have the form GiR-{()'^G2 with k > 1 and some operators 
Gi,G2. These terms vanish after integration along Fq. 

Let us consider the structure of all remaining operators of order lower 
than —N again. They will contain b and its derivatives as factors. The 
only term where only one factor b is present, as it follows from (4.3), 
will be the one of the form 2b(P_ — C)~^! obtained by commuting the 
first and second factors in (P_ — C)(2b)(P„ — (). The integral of this 
term along the contour vanishes. All the remaining terms will contain 
at least two entries of b or at least one derivative of b. By Proposition 
2.1, used with proper /, this means that each such term is bounded as 
acting into the weighted space with weight {x)^^~^^. This boundedness 
property is preserved after the contour integration. 

The detailed combinatorics in the above reasoning as well as the 
inductive procedure enabling us to pass from ^ = 1 to an arbitrary q, 
is explained in [15]. □ 

Thus the operator &g is close to the projection Pg. Together with 
the closeness of (5q to Qq, this reasoning proves Theorem 4.1. 

Now we add a perturbation by a smooth electric potential V{x) G Sp. 
Similar to [15] under such perturbation the spectral subspaces 'almost' 
do not change. Note, first of all, that the perturbation of P_ by V 
is relatively compact, therefore, again, the spectrum of the operator 
P_{V) = P_ + V between Landau levels is discrete. We can change 
the contours Tg a little, so that they do not pass through the eigenvalues 
of P_(y). We denote by the spectral subspaces of P-(T^) corre- 
sponding to the spectrum inside Tg, by the corresponding spectral 
projections, and by RY^{C) the resolvent of P_(V^). 

Proposition 4.5. The projections P^ and Qg are close in the sense 
used in Theorem 4-1- Moreover, {x)^^^^T{P^ — Qq)T' is hounded for 
any finite products T, T' of creation and annihilation operators. 

Proof. We will prove that the projection Pg is close to Pg, then the 
result will follow from Theorem 4.1. 

We use the representation of projections Pg and Qg by means of 
resolvents and subtract: 



{x)-^^'P^ -Pg^ -{2ni)-\x)-^^' / R^OVR'^iOdC = (4.6) 
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The last term in (4.6) is bounded and stays bounded after the mul- 
tiphcation by the creation and annihilation operators, as soon as N 
is large enough, by Lemma 4.2. With the leading terms in (4.6), we 
can perform the same procedure as when proving Proposition 4.4. We 
commute the resolvent with V and with the terms arising by com- 
mutation and so on, aiming to collect the resolvents together all the 
time. We arrive at a number of terms of sufficiently negative order, 
thus bounded before the integration, and terms with all entries of the 
resolvent collected together, and thus vanishing after the integration. 
The boundedness property with weight follows again, as in the proof of 
Proposition 4.4, from the fact that after the commutation, all surviving 
terms will contain at least two entries of V or an entry of a derivative 
of V^. □ 

5. Spectrum of Toeplitz-type operators 

We move on to the study of the sphtting of Landau levels of our 
operators. Similar to [13], [10], [15], [7], the properties of the eigenval- 
ues of the perturbed operators arc determined by the properties of the 
spectrum of certain Tocplitz-typc operators. 

Usually, by Toeplitz operator one understands an operator of the 
form %{W) — PWP where P is the orthogonal projection onto some 
subspace Q in L2 and W is the operator of multiplication by some 
function. Alternatively, if we consider the Toeplitz operator as acting 
in Q it can be written as T(iy) = PW. Usually, the subspace Q 
consists of functions, related with analytical ones, for example. Hardy 
or Bergman spaces. 

In [13], [10] such operators were considered, with Q being one of 
Landau subspaccs and the electric potential V acting as W. We will 
study, as in [15], this latter kind of Toeplitz type operators, with the 
space Tio of zero modes of the perturbed Pauli operator acting as Q 
and some differential operator acting as W. The result below differs 
essentially from the one in [15]. 

Let V, b, be real functions in S13, (3 < —2. We consider the Toeplitz- 
type operator in Ho' 

%ou = %o{y)u = ^oQ'(P- - Ag + V^)Q%, u e Ho- (5.1) 
This operator, by (2.10), corresponds to the quadratic form 

t^[«] = ((P_ - A, + \/)Q^u, Q%) = (P+ - A,+i)IIQ'+'w|l' 

-Ag+i||Q%f + {{V - 2b)Q^M, Q%) ,ueno. (5.2) 

We are going to study the spectrum of To(V^) . We denote by 
— A^(To(V^)) positive, resp., negative eigenvalues of the oper- 
ator %o{V). The distribution functions n±(A), A > 0, are defined 
as n±{X) — n±(X;%o(V)) — #{n- '■ > A}. The s-numbers 

Sn of the operator %o{V) are just the absolute values of A J ordered 
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non-increasingly, and their distribution function equals n{X,%o{V)) = 
n+{X,%o{V)) + n4X,%o{V)). 

We denote by V = V^fV, b] the effective weight 

V[V, b] = C,(V + 2qh), C, = q\(2By, (5.3) 

and set for any function W E±{X,W) = (27r)-^B° meas {x e : 
±W{x) > A}. 

Proposition 5.1. For the eigenvalue distribution function n±{X; %q{V)) 
of the operator 'Zo{V) the following estimate holds: 

n±(A;To(y)) < CaI, A^O (5.4) 

Moreover, suppose that the function E±{X,Y) (for one or for both 
signs) satisfies the estimate 

E±{X,Y)>C'X^,X^O (5.5) 
and the regularity condition 

hmhmsup ^ . = 1. (o.d) 

Then for the distribution function n±{X,%Q{y)) of the eigenvalues of 
To(^) the asymptotic formula holds 

n±(A,2:o(V^))~i?±(A,V), A^O. (5.7) 

The conditions of the form (5.5) and (5.6) are traditional in the study 
of asymptotic properties of operators with (possibly) non-power behav- 
ior of eigenvalues, see, e.g., [11], [7], [8]. The first of them indicates that 
the 'rate of decay' (5 is chosen sharply, so that various remainder terms 
are, in fact, weaker than the leading one. The second condition means 
that the function i?±(A, V) grows sufficiently regularly. It enables the 
use of various kinds of perturbation techniques. 

We note also that for = and V — the effective weight V 
vanishes. This property corresponds to the fact that P_ has TYo as its 
space of zero modes. 

Proof. By Proposition 3.4, the expression (5.2) can be transformed to 

^[u] — (Vii, u) -\- (vu, u), u eHo, 

where v G S/s^s and V = Vg[V^, b] e Sp. By Lemma 2.5 in [7], for the 
operator defined by this quadratic form, the asymptotic relation 

n±{X) ~ (27r)-^B° meas {x e : V(a;) + v{x) > A}, A ^ 0, 

holds, as soon as the measure on the right-hand side grows not slower 

2 

than A/3 and is regular, i.e., if the conditions (5.5), (5.6) are satisfied. 
This reasoning takes care of the second part of the Proposition. The 
first part follows from Proposition 2.3 in [7], without qualified lower 
estimate. The function v, due to its decay, does not contribute to the 
main order of the eigenvalue estimates and asymptotics. □ 
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Now we can establish the spectral estimates and asymptotics for a 
similar Toeplitz operator on the spectral subspace corresponding to the 
cluster around the Landau level for an arbitrary q. 

Proposition 5.2. Let b,V E Sp, P < —2. Consider the Toeplitz-type 
operator 

%{V) = P,(P_ - A, + V)P,. (5.8) 

Then for the distribution function of s-numbers of the eigenvalues of 
^qiy) the estimate holds 

n(A,T,(\/)) = 0(Al). (5.9) 

//, moreover, the effective weight V defined in (5.3) satisfies the con- 
ditions (5.5), (5.6) then the asymptotic formula 

n±(A, %,{V)) - £;±(A, V + 2qh) (5.10) 

holds. 

Proof. By Proposition 4.4, the spectral projection Pq can be approxi- 
mated by &q = C^Q^PoQ^- We use (4.2) to obtain 

%,{v) = e,(p_ - A, + v)eq + (p,(p_ - A,)p, + zj(p_ - A, + v)e, 

+6,(P_ - A, + y)(P,(P- - A,)P, + Z,)+ (5.11) 
(P,(P_ - A,)P, + Z,)(P_ - A, + y)(P,(P_ - A,)P, + Z,). 
We re-arrange (5.11) to separate the leading terms: 

%,{v) = e,(p_ - A, + T/)e, + p,(p_ - A, + i/)p,(P- - A, + y)p,e5 
+egP,(P- - Ag + y)p,(P- - Ag + y)p,+ 

P,(P_ - A, + y)P,(P_ - A, + 1/)P,(P- - A, + V)P, + Y. 
So, we have 

%,{v) = 6,(p_ - A, + v)&, + T,(y)2e, + 6,x,(i^)2 + %^{yY + y. 

(5.12) 

Let us consider the structure of the remainder operator Y. On the 
one hand, it contains terms having the factor Zg. By Proposition 4.4, 
the operator (x)^'^+*(P+)"^Zg is bounded for N that can be chosen 
arbitrarily large. Thus, by Proposition 2.2, the s-numbers of Z^ satisfy 

2 

the estimate n(A, Z^) = o{Xi^), A ^ 0. Therefore all terms in Y 
containing Z^ satisfy the same kind of s-numbcrs estimate. 

The terms in Y not containing Z^, contain as factors the function 
V , the projection Pq and some compact operators that remain com- 
pact after the multiplication by any power of P+. For the operator 
y(P+)-^, by Proposition 2.2, the estimate n(A,y(P+)-^) = 0(A?) 
holds. After the multiplication by a compact operator, the O symbol 
can be replaced by a in this formula. As a result, we obtain 

n(A,Y) = o(A^),A^O. (5.13) 
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Now consider the second, third, and fourth terms on the right-hand 
side in (5.12). They contain the square of the operator %q{V) multiphed 
by some bounded operators. Therefore the rate of decay of s-numbers 
of these terms is faster than the rate of decay of the ones of %q{V). 
So, the estimate (5.9) or the asymptotic formula (5.10) will be proved 
as soon as we estabhsh these formulas for the operator D = 6g(P_ — 
Aq + V)&q. The quadratic form of this operator is (Dm, m) = ((P- — 
Aq + V)&qU, &qu). Recalling the definition of &q in (4.1), we have 

(Du, u) = Cg(Q«PoQ'^i, (P- - A, + V)Q'iPoQ''u). 
We set here v — PqQ'^u, v e Ho, to get 

(Du, u) = C^'iQ'^v, (P_ - A, + V)Wv), V e Ho. (5.14) 

The spectral estimate and, under the conditions (5.5), (5.6), the 
asymptotics for the operator defined by the latter quadratic form is 
given by Proposition 5.1, where the factor C^^ is responsible for the 
replacement of V by C^^\ = V + 2qh. This result carries over to the 
quadratic form (Du, u) using the fact that u = Cg(l + K)Cl'^v for some 
compact operator K, which is explained in [15]. □ 

The operators of the above type will be used in the next Section in 
order to find the leading term in the eigenvalue asymptotics in clusters. 
Another type of Toeplitz operators will be needed in order to perform 
a block digitalization of the Pauli operator. 

Proposition 5.3. Let W be a function in Sj3,(3 < —2. Consider the 
operator %'{W) = (1 — Pq)WPq. Then for the distribution function of 
singular numbers of%'(W) 

n(A, %'{W)) = o(At), A ^ 0. (5.15) 

Proof. The quadratic form of the operator %' (WyX' (W) equals 

(l\W)u, %'{W)u) = ((1 - Pq)Wu, (1 - Pq)Wu) (5.16) 

= {[W, Pq]PqU, [W, Pq]PqU),U G Hq. 

Consider the case g = first. By Proposition 2.2 the operator \W, Pq\ 
can be represented as 

[W,Po]=L{xf-' (5.17) 
for some bounded operator L. By Proposition 5.1, the operator \W, Pq] 
has singular numbers with the required decay rate. 
For an arbitrary > by Proposition 2.1, 

[W,Pq]^L{x)P-'; (5.18) 

We use the approximation of the projection Pq found in Section 4: 

[W, Pq]Pq = L{xY-'eq + L{xf-\Pq - S,), (5.19) 

where &q is defined in (4.1). The spectral estimate for the first, leading 
term in (5.19) involves the projection Pq and the task of estimating its 
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singular numbers reduces to the already considered case q = 0, see 
(5.14). By Proposition 4.4, the second term in (5.19) can be written as 
L{x)f^~^P^^ K with as large N as needed and a compact operator K. 
The required spectral estimate for this operator follows now from the 
second part of Proposition 2.2. □ 

6. Perturbed eigenvalues 

Now we are able to establish our main result about eigenvalue asymp- 
totics and estimates for the perturbed Schrodinger and Pauli operators. 
For a sclf-adjoint operator L, wc denote by N{X,fi) = N[X,iJ,]L) the 
number of eigenvalues of L in the interval (A,/x). 

Theorem 6.1. Let V, b e Sp,(3 < —2. Fix an integer q > and let 
X± be some fixed real numbers, \± ^ Ag, |A± — Ag| < B°, such that 
they are not the points of the spectrum ofP- + V . 
For P_(y) the following estimates hold for A — > 0+ 

N{K, + A, A+; P_(y)) = O(A^), 7V(A_, A, - A; P_(y)) = 0(A§). 

(6.1) 

// the effective weight Vq[y, b] = CqiV + 2gb) or, what is equivalent, 

V + 2qh satisfies the conditions (5.5), (5.6) with the sign "+" then, 
asymptotically, 

N{Ag + A, A+; P-{V)) ~ E+(A, V + 2qh), A ^ 0, A > 0. (6.2) 

If the effective weight satisfies the conditions (5.5), (5.6) with the sign 
then, asymptotically, 

A^(A_,Ag-A, ;P_(y))~£;_(A,y + 2gb), A ^ 0, A > 0. (6.3) 

Similar results hold for the operator 'P j^iV) and for the Schrodinger 
operator H(K) with the following obvious modifications, corresponding 
to (3.3). For the Pauli operator P+(\/) one should replace in the es- 
timates of the form (6.1) and in the asymptotic relations of the form 
Aq by Ag+i, A± by X± + 2B° and V by V + 2b. For the Schrodinger 
operator il{V) one should replace Aq by Aq + B°, A± by X± + B° and 

V byV + h. 

Proof. We will use the following statement about the eigenvalue distri- 
bution of perturbed operators. If Lq, Li are two self-adjoint operators, 

moreover, Li is compact, then 

N{ni, Lq + Li) < N{ij,i - Ti, + T2; Lq) + n(ri; Li) + n(r2, Li) 

(6.4) 

for any interval (/Xi,/i2) and any positive numbers Ti,T2. The proof of 
(6.4) can be found, e.g., in [11], Lemma 5.4. 

We show now the upper asymptotic estimate in (6.2). For a fixed q, 
we take as Lq the operator 

Lo = Pq(P- + V)Pq + (1 - P,)(P_ +V){1- Pq), (6.5) 
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and as Li the operator 

Li^{l-P,)VPg + P,V{l-P,), 

so that Lq + Li = P_ + V. We fix some e > 0, and apply (6.4) for 
/ii = Ag + A, A+, Ti = r2 = eA. The spectrum of Lq is the union of the 
spectra of the summands in (6.5). The asymptotics of the eigenvalues 
of Pq(P- + V)Pq is given by the Proposition 5.2. The second term 
in (6.5) contributes to the spectrum near only with finitely many 
points. Thus, iV(Ag + -ri, A^ + /i2 + rg; Lq) ~ E+{\{1 - e);V + 2qb). 
On the other hand, for the spectrum of the operator Li, by Proposition 

2 

5.3, we have n{e\,Li) = o{{eX)i^). We substitute these asymptotic 
estimates into (6.4), divide by E+{X,V + 2qh) and pass to limsup as 
A — > 0. We arrive at 

limsup ^ , , , — - — „ , , < 1. 

A-.0 E+{X{l-e);V + 2qh) ' 

Due to the arbitrariness of e, by our assumptions, this implies the upper 
asymptotic estimate in (6.2). 

All other upper estimates in the Theorem are established analo- 
gously. The lower asymptotic estimate in (6.2), (6.3) is estabhshed 
in the same way, just interchanging Lq and Li in (6.4). □ 
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